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Abstract 

The ordering matrix, which was originally introduced by de Gennes, is a well-known 
mathematical device for describing orientational order of biaxial nematic liquid crys- 
tal. In this paper we propose a new interpretation of the ordering matrix. We slightly 
modify the definition of the ordering matrix and call it the geometric order parameter. 
The geometric order parameter is a linear transformation which transforms a tensorial 
quantity of an individual molecule to a tensorial quantity observed at a macroscopic 
scale. The degree of order is defined as the singular value of the geometric order 
parameter. We introduce the anisotropy diagram, which is useful for classification 
and comparison of various tensorial quantities. As indices for evaluating anisotropics 
of tensorial quantities, we define the degree of anisotropy and the degree of biaxiality. 
We prove that a simple diagrammatic relation holds between a microscopic tensor and 
a macroscopic tensor. We provide a prescription to formulate the Landau-de Gennes 
free energy of a system whose constituent molecules have an arbitrary shape. We 
apply our prescription to a system which consists of D2/t-symmetric molecules. 
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1 Introduction 



A system which consists of asymmetric molecules can exhibit various noticeable phenomena. In par- 
ticular, biaxial nematic liquid crystals are interesting subjects of current research. Biaxiality means 
that physical properties of an individual molecule or an ensemble of molecules are not invariant under 
any rotations. On the other hand, uniaxiality means that properties of a molecule or an ensemble 
of molecules are invariant under rotations about a fixed axis. Anisotropy is a general concept which 
implies either uniaxiality or biaxiality. Although most of molecules are not uniaxially symmetric in 
a rigorous sense, it is possible that some properties of molecules are effectively uniaxial. In a simple 
nematic liquid crystal uniaxial molecules are aligned to exhibit a uniaxial order at a macroscopic scale. 
However, in a more complex system under some circumstance, it may happen that uniaxial molecules 
exhibit a biaxial order, or it may also happen that biaxial molecules exhibit a biaxial order. 

Biaxial nematic phases have been intensively studied in liquid crystal physics. Williams pQ noticed 
biaxial anisotropy in optical properties of nematic liquid crystals in a magnetic field. Taylor, Fergason, 
and Arora [2] found a biaxial smectic C phase. Freiser [3] began a theoretical study of phase structures 
of asymmetric molecules and predicted the existence of a biaxial nematic phase. Alben [1] calculated 
the Landau free energy of a system of biaxial molecules and predicted the existence of a biaxial 
phase. Straley [6] introduced four order parameters, (S, T, U, V) in his notation, to describe nematic 
order structures of an ensemble of molecules which have the point group T>2h as their symmetry, de 
Gennes and Prost 7J introduced a set of generalized order parameters, which is called the ordering 
matrix, to describe nematic order structures of molecules which have an arbitrary shape. Yu and 
Saupe 8J observed a biaxial phase in an experiment and obtained a phase diagram in the concentration- 
temperature coordinate. There the biaxial phase appeared between two distinct uniaxial phases. 
Boonbrahm and Saupe [§] studied effects of temperature and magnetic field on a thin film of biaxial 
nematic liquid. Allender, Lee, and Hafiz [Till E] constructed the Landau-de Gennes free energy of 
Z^-symmetric molecules up to the sixth order in the Straley variables (S, T, U, V). Bunning, Crellin, 
and Faber [13] studied experimentally an effect of molecular biaxiality on bulk properties, particularly 
on the magnetic anisotropy. Gramsbergen, Longa, and de Jeu [14j wrote a review on the Landau 
theory of the nematic-isotropic phase transition. However, the effects of biaxiality of molecules were 
not sufficiently considered in their review. Remler and Haymet [T5] gave a complete formulation for 
describing interactions between asymmetric molecules and applied their formulation to the analysis of 
the Landau free energy. Mulder [TO] formulated a model of sphcro-platelet molecules which interact 
by exclusive volume effect. Since his model has the T>ih symmetry, the order is described by the four 
order parameters of Straley. Solving the problem by the mean field approximation, he showed that a 
transition between the isotropic phase and the biaxial phase can occur. 

The discoveries of a thermotropic biaxial phase in a system of bent molecules (banana-shaped or 
boomerang-shaped molecules) by Madsen, Dingemans, Nakata, and Samulski [5U] and by Acharya, 
Primak, and Kumar [5T] renewed the interest in biaxial nematics [55]. In their experiments it was 
observed that biaxial molecules exhibited biaxial orders without application of external fields nor 
boundary effects. Their discoveries have been stimulating intensive researches in this field [23 | 154 ] [251 , 
26, 28, 29] . Merkel et al. [53J measured biaxiality parameters by infrared absorbance measurements and 
compared the observed data with a result of the Landau-de Gennes model. Bates and Luckhurst [53] 
studied the phase diagram of a liquid which consists of V-shaped molecules by the Monte Carlo 
simulation. They showed existence of biaxial phases in the diagrams whose coordinates are temperature 
and various anisotropy parameters like the bending angle of the molecule. 

However, for a theoretical analysis of biaxial nematic phases, it seems that there is still a confusion 
in descriptions of anisotropies. In other words, it is necessary to invent a more useful and comprehensive 
method for describing anisotropies of molecules and nematic phases. Let us discuss issues which exist in 
the present method for describing anisotropies of general nematics. For a liquid crystal which consists 
of uniaxial molecules, a well-known device to characterize orientational order of nematic phase is the 
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tensorial order parameter 

A = (n®n—'^iy (1) 

Here n is a unit vector which is fixed along the axis of each molecule and viewed from a laboratory 
observer. Since the molecules execute thermal motion, the direction of the vector n fluctuates. The 
brackets (• • •) mean a statistical average. The components of the tensor A are written as 



Aij = (mrij - -Sij), i,j = 1,2,3. (2) 



The eigenvectors and eigenvalues of A indicate alignment of the molecules. The matrix A can be 
diagonalized and parameterized as 



-a + r 
-a-T 
2a j 



A=\\ II -n - t | . (31 



When a = r = 0, the system is in an isotropic phase. When a ^ and t = 0, the system is in a 
uniaxial nematic phase. The value of a is in the range — \ < a < 1. When t ^ 0, the system is in a 
biaxial nematic phase. 

If the molecule itself is biaxial, we may introduce another order parameter, 



B = (l®l-m®m). (4) 

Here I, m, n are mutually orthogonal unit vectors fixed on each molecule. The quantity B characterizes 
the biaxial anisotropy of the nematic phase. In general, we may associate principal values Ai,A2,A3 
with the principal axes I, m, n and define the order parameter 

C = (^Xil ® I + A 2 m ® m + \ 3 n ® n - ~(Ai + A 2 + A 3 )i^. (5) 

These order parameters, A, B, and C, are useful for characterizing orientational order structures 
of nematic phases. But there are several difficulties in their application to molecules which have 
an arbitrary shape. First, there is no a priori reason to choose the molecular axes l,m,n for an 
asymmetric molecule. If the molecule is rectangular, choice of the axes is rather obvious. However, for 
a molecule which has no symmetry, choice of the axes is not unique. There are various candidates for 
the molecular axes; we may take the principal axes of the inertia tensor, the dielectric susceptibility 
tensor, the electric quadrupole tensor, or the magnetic susceptibility tensor of the molecule. In general, 
the axes defined by them do not coincide. Thus, there is no unique definition of the molecular axes. 
Second, distinction between the uniaxiality and the biaxiality becomes ambiguous since the eigenvalues 
of A and B depend on the choice of the molecular axes. Furthermore, there is no reason to choose 
a unique set of the principal values Ai, A2, A3 in the definition of the tensor C. Third, the relation 
between the anisotropy of a molecule and the anisotropy of a macroscopic phase is vague in this kind of 
analysis. It can happen that uniaxial molecules exhibit a biaxial phase. It is also possible that biaxial 
molecules exhibit a uniaxial phase. Thus a systematic method to compare the molecular anisotropy 
and the macroscopic anisotropy is desirable. 

de Gennes [7] introduced the ordering matrix 

Sabij — I (^Rai Rbj ~ -^SabStj^ , a, 6, i, j = 1, 2, 3 (6) 

for characterizing alignment of molecules in a nematic phase. Here R a i = L a ■ Mi is an inner product 
of the laboratory orthogonal frame (Li, L 2 , L3) with the molecular orthogonal frame (Mi, M 2 , M3). 
The symmetrized tensor 



Gabij — { — (Rai Rbj + Rbi Raj) — g<5afc Sij j (7) 
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is more useful and meaningful as will be shown in this paper. Although the ordering matrix is applicable 
to molecules of an arbitrary shape, it is still difficult to read out geometrical and physical implications 
from the ordering matrix. 

In this paper we introduce a new approach for characterization and analysis of anisotropics of a 
molecule and a bulk phase. However, here we describe the outline of this paper. In our discussion, the 
adjective microscopic means intrinsic properties or quantities which an individual molecule possesses. 
On the other hand, macroscopic means average properties or quantities observed in an ensemble of a 
large number of molecules. If each molecule has a tensorial quantity Uj and if the molecule changes 
its direction, the tensor is transformed to t ao — ■ R a i Rbj Uj by a rotation matrix R a i. We assume 
that Uj is a traceless symmetric tensor. The quantity observable at a macroscopic scale is a statistical 
average (t a b) = J2i j(Rai Rbj) Uj = X)i jG a UjUj- This is an equation defining the geometric order 
parameter G a bij- Thus, the geometric order parameter can be regarded as a bridge which relates 
the microscopic quantity Uj to the macroscopic quantity (t a b). Since the geometric order parameter 
G = (G a bij) is a linear transformation t i— > Gt, its property is completely analyzed by the method of 
singular value decomposition. In Sect. [2] we will introduce the geometric order parameter and discuss 
its properties. 

After understanding the geometric order parameter, the remaining task is to characterize anisotropics 
implied by the individual tensors, Uj an d (tab)- To visualize the anisotropic property of a tensor we 
introduce an anisotropy diagram, in which each tensor is represented as a point in a plane. Then we 
define the degree of anisotropy a(t) and the degree of biaxiality j3(t) of the tensor t — (Uj). Sect. [3]is 
an introductory discussion for providing the indices of anisotropy and Sect. 2] is an explanation of the 
anisotropy diagram. 

Furthermore, the geometric order parameter enables us to compare anisotropies of the microscopic 
tensor t and the macroscopic tensor (t) = Gt. We found that in the anisotropic diagram there is a 
simple geometric relation between the microscopic tensor and the corresponding macroscopic tensor. 
In Sect. Owe will prove some theorems to ensure the micro- macro relation. This section is a highlight 
of this paper. 

In Sect. [H] we will show simple applications of our method. In Sect. [7] we restrict our consideration 
to molecules which have the D^h symmetry. Then, we will reproduce the four order parameters of 
Straley. In Sect. [5] we will give a general prescription to formulate the Landau-de Gennes free energy 
model. There we refer to the theorem which tells a complete set of ingredients of the Landau-de 
Gennes free energy. In the appendix we prove the theorem. A real molecule may have various tensorial 
quantities which are not simultaneously diagonalizable. Our prescription is applicable even to such 
a general system. Finally, we apply our prescription and obtain a complete Landau-de Gennes free 
energy for the r>2/i-symmetric molecules. Sect. |9] is devoted to concluding remarks. 

We would like to emphasize that our method for characterizing anisotropies is applicable to a general 
system in which molecules may have arbitrary shapes and arbitrary tensorial quantities. Our method 
is systematic and unambiguous. The anisotropy diagram will help both qualitative and quantitative 
understandings of anisotropies. Our prescription for formulating the Landau-de Gennes free energy 
enables us to construct a complete invariant polynomial which contains neither too many nor too few 
terms. 

2 Geometric order parameter 

In this section we introduce the geometric order parameter. Although it is just a modified version of 
de Gennes' ordering matrix, it will give a clear and new interpretation of the ordering matrix. 

Assume that a molecule has an intrinsic vectorial quantity v = (i>i, ua, U3), which can be, for 
example, an electric dipole moment. When the molecule rotates, the vector v is transformed to 

3 

v = Rv, or v a = 22 Rai Vi (8) 

»=i 
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by a three-dimensional orthogonal matrix R G 5*0(3). The matrix elements R a t satisfy ^ a R a iR a j = 
5ij and ^2nR a iRu = $ab- Each molecule can be transformed by a different rotation matrix. Since 
liquid crystal is an ensemble of molecules, the quantity observed in the laboratory is the average 



3 

E 

i=l 



{Rai)Vi. 



(9) 



Once we know the matrix elements (R a i), we can calculate the average (w) = (R)w for any vectorial 
quantity w of the molecule. Most of liquid crystals have no polarity and hence (R a i) are usually zero. 

Next, assume that the molecule has an intrinsic tensorial quantity t = (£y), which may be a 
dielectric susceptibility or an electric quadrupole moment. When the molecule rotates, the tensor t is 
transformed to 



i=(R® R)t, or 



tab 



(10) 



^ Rai Rbj t'ij • 

Any tensor t can be decomposed into the scalar component, the antisymmetric component, and the 
traceless symmetric component as 



3 



6 k=l 



tkk 



-(Uj + tji) - -5ij 

fc=l 



(11) 



If the tensor t is traceless and symmetric, the transformed tensor i is also traceless and symmetric. 
Hence we can write the components of t as 



1 - 1 3 3 

tab = y(tab + ha) ~ o<5ab t C c = 



-{Rai Rbj + -Rbi ^aj) — g^ab 8ij 



Thus, the transformation law of traceless symmetric tensors is described as 



(12) 



tab ^ Qabij tij 



with the symmetrized traceless matrix 



Qabij — 2 (Rat Rbj ~t~ ^fei Raj ) 3^°^ * 

Then the average, which is an observable at a macroscopic scale, is given by 



(13) 



(14) 



(t) = (Q)t = Gt, Or (i„i,) = (Qabij)tij = J2 Ga 



bij tij 



(15) 



The defining equation of (Q a bij) — Gabij is Eq. ([7]). Once we know the matrix elements (Qabij), we 
can calculate the average (u) — (Q)u for any tensorial quantity u of the molecule. It is not necessary 
that the tensors t and u have common principal axes. We call (R) and (Q) the geometric order 
parameters. More specifically, we may call (R) the geometric order parameter for vectors while we call 
(Q) the geometric order parameter for traceless symmetric tensors. In our approach, the macroscopic 
observable (t) is calculated as a product (t) = (Q)t of the geometric order parameter (Q) with the 
molecular intrinsic quantity t. In this treatment we can analyze anisotropies of t and (t) separately. 

Superficially the ordering matrix G a uj = {Qabij) has 3 4 = 81 components but actually it has only 
25 independent components [5]. The geometric order parameter G transforms a traceless symmetric 
tensor t into another traceless symmetric tensor (t) = Gt. The set of all traceless symmetric tensors 
forms a 5-dimensional vector space and G is a linear transformation of the space of traceless symmetric 
tensors. Hence the ordering matrix G has 5 2 = 25 independent components. This fact can be verified 
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also by counting independent components of G a bij which are restricted by the traceless and symmetry 
conditions 

3 3 

^ ' Gggij — 0, y ' Ggbii — 0, G a bij — Gbaij — G a bji- (16) 

a—1 i—1 

We would like to have a representation of the geometric order parameter in which only independent 
components appear explicitly. For this purpose we will introduce the reduced ordering matrix G^ v in 
the following. First, we define an inner product of two tensors t and u as 

3 

(t, u) = Tr{t T u) = ^ U, (17) 
Here t T is the transposition of t. It is allowed to make a product tu of two tensors as 

3 

(tu) ik = }JijUj k . (18) 
j'=i 

Second, we introduce a basis • • • ,£ 5 } of the space of traceless symmetric tensors, 

^0 1 0\ 1 /0 1\ 1 /° ° N 

€ x = 4= I 1 , £ 2 = -= , £ 3 = -= 1 



\0 0/ ~ 0/' x - 

'1 0\ , /-l ON 



4 1 0, 



C * = "75 -1 ' Cs = IS - 1 I ■ ( 19 ) 
V2 \0 0/ V6 V 2y 

They satisfy (£„,£„) = i^i, with respect to the inner product (jTTJ) . We write the components of £ M 
as £/i y with indices fi — 1, • • • ,5 and i,j = 1, 2, 3. An arbitrary traceless symmetric tensor t can be 
expressed as a linear combination of • • • , £ 5 }, 

5 

* = E c ^m ( 20 ) 

with the coefficients c M = (£ M ,t). Finally, we define a 5-dimensional matrix G = (G M „) by 

5 

G^ = ^^G^, i/=l,.-.,5. (21) 

These matrix elements G MI , can be calculated as 

G^ = G£„) = ^ £/j,,ab G a bij £,v,ij- (22) 
a, 6, i, j — 1,2,3 

We call the 5-dimensional matrix G = (G^j,) the reduced ordering matrix. By the definition it has 
5 2 = 25 independent components. When the macroscopic tensor (t) = Gt is expanded in the basis as 
(t) = 53u=i ^m^m' ^ s com P onen ts are given by 

5 

-i^^^c. (23) 
i/=i 

The components G a uj of the original geometric order parameter can be reconstructed from the com- 
ponents G M „ of the reduced ordering matrix as 

5 

Gabij ^ y £,[i,ab G ^ u £,v,ij • (^^) 
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Therefore the reduced ordering matrix contains the same information as the geometric order parameter. 

To read out the implication of the geometric order parameter we apply the singular value decom- 
position on it. Here we review the definition of the singular value decomposition of a matrix. For a 
matrix G = (G^) if a set of vectors c a — {c^ a ), d a — (e? MQ ) and real numbers a a (/x, is, a = 1, • • • , 5) 
satisfy 

5 5 5 5 

u— 1 /j— 1 1 M— 1 

then the vector c a is called the right singular vector, d a is called the left singular vector, and a a is 
called the singular value. The above equations can be written more concisely as 

Gc a = a a d a , G T d a = cr a c a , (c a ,cg) = S a/3 , {d a ,dp) = 5 a p. (26) 

Here G T is the transposed matrix of G. It is always possible to make o~ a non-negative by choosing c a 
and d a suitably. For a symmetric matrix G — G T , the left singular vector and the right singular vector 
coincide and they are called an eigenvector. In this case the singular value is called an eigenvalue. From 
the singular vectors and singular values we can construct matrices 

C = (c pQ ), D = (dpa), S = (a a 5 a/ 3). (27) 

Note that S is a diagonal matrix. Then the set of equations (|25p is equivalent to 

GC = L>£, D T G = SC T , C T C = I, D T D = /, (28) 

which implies D 7 GC = E. This is a generalization of diagonalization of a matrix. It can be rewritten 
as 

G = DYiC T (29) 

and this expression is called the singular value decomposition of G. 

Now we apply the singular value decomposition to the reduced ordering matrix G to understand 
the implication of the geometric order parameter. Once we know the singular vectors of G, c a — (c^ a ) 
and d a = (d Ma ), we can construct tensors 

5 

Then the definition of singular vectors ((2"6f implies 

Gt a =a a u a , (t a ,t/3) = 5 a p, (u a ,up) = Sap. (31) 

On the other hand, as discussed at (fT5|) . when the molecule has an intrinsic physical quantity t a , 
the average (t a ) = Gt a will be observed by a macroscopic measurement. The observed value is now 
given as (t a ) = a a u a . The coefficient a a takes its value in the range < a a < 1 and is called the 
degree of order or the strength of realization. The reason why o~ a is in the range < o~ a < 1 will 
be explained in Sect. [5] as a corollary of theorem 1. The tensor t a is called the microscopic singular 
tensor and u a is called the macroscopic singular tensor. It is convenient to arrange them in the order 
ci > C2 > • • • > (T5. Then, if each molecule has a quantity represented by t a , the ensemble of molecules 
exhibits the quantity u a at the macroscopic scale with the strength a a . If a a = 0, the effect of the 
molecular quantity t a disappears at the macroscopic scale. 

Let us summarize the above discussion. The equation (|15j) relates the microscopic tensorial quantity 
t to the macroscopic observable (t). The equation (t) = Gt can be rewritten symbolically as 

(macroscopic observable) = (geometric order parameter) x (microscopic quantity). (32) 



u, 



u=l 



fidp, a . 



(30) 
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Furthermore, the equation (|31[) tells that the molecular quantity t a manifests itself as the macroscopic 
quantity u a with the strength a a . This relation Gt a = a a u a can be expressed symbolically as 

(geometric order parameter) x (microscopic singular tensor) 

= (strength of realization) x (macroscopic singular tensor). (33) 

In this way we can read the implication of the geometric order parameter G. 

We would like to mention another interesting property of the geometric parameters. In a nematic 
phase orientations of molecules are fluctuating. The orientation of each molecule is specified with 
a three-dimensional rotation matrix R £ 50(3). Then distribution of the molecular orientations is 
described by a probability density function p(R) over 50(3) and the average of a physical quantity 
f(R) which depends on the orientation of a molecule is given by the integral 

(/) = / f(R) P (R)dR = ^J f(R)p(R) sin0<20##- (34) 

In the last line we used the Euler angles (6, <fi, -0) to specify the rotation matrix R. Note that Q a Hj (R) 
defined in (|14|) is a function of R E SO (3). Furthermore, if we define 

,ab Qabij 7 (35) 

a, j— 1,2,3 

Q\tv{R) is also a function of R G 50(3). The matrix Q(R) = {Q^{R)) forms a 5-dimensional 
irreducible representation of the rotation group 5*0(3). Namely, it satisfies Q(RR') = Q(R)Q(R') for 
any R, R' 6 50(3). If we know the probability density p(R), we can calculate the averages (R a i) and 
(Qfiu). Actually, the inverse of this statement holds. Once we know the averages (R a i) and (Q^), we 
can determine the probability density p{R) via 

3 5 

p(R) = 1 + 3^ (R ai ) R ai + 5 @i>») + ■■■ ■ ( 36 ) 

a,i— 1 

This equation is regarded as an expansion of p(R) in powers of R a i- It is easily proved by applying the 
Peter- Weyl theorem |19j of group representation theory. In this way the geometric order parameters 
completely characterize the geometric and statistical properties of the nematic phase. 



3 Elementary attempts to characterize anisotropy 

In the previous section we argued that the microscopic tensorial quantity t is related to the macroscopic 
tensorial quantity (t) via the geometric order parameter G as (t) — Gt. We also showed that the 
implication of the geometric order parameter G can be analyzed via the singular value decomposition. 
The remaining problem is to provide a systematic method to analyze properties of tensorial quantities, 
t or (t), particularly their anisotropy. This is the subject we will discuss in this section. 

Here we discuss briefly some attempts to characterize anisotropy of a symmetric tensor t — (tij) 
(the trace is not necessarily zero). The tensor has three principal axes and three eigenvalues Xx, A2, A3. 
By choosing the spacial coordinate suitably, we can transform it in a diagonal form 

/A a \ 

t = A 2 . (37) 

\o A 3 y 

When the three eigenvalues coincide, it is said that the tensor is isotropic. When two of the three 
eigenvalues coincide, the tensor is uniaxial. When the three are distinct, the tensor is biaxial. We 
would like to define indices which indicate quantitatively the degree of anisotropy and the degree of 
biaxiality. 
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As a candidate for the index of anisotropy we may introduce 

A = (Ai - A 2 ) 2 + (A 2 - A 3 ) 2 + (A3 - A1) 2 . (38) 

It is obvious that A is non-negative. If and only if A = 0, the tensor is isotropic. On the other hand, 
we define the average of the eigenvalues m = ^(Ai + A2 + A3) and the standard deviation 

A = 1 1 (Ax - m) 2 + (A 2 - m) 2 + (A3 - m) 2 }. (39) 

We may take A as another index of anisotropy but actually they are related as 

A = 9A. (40) 

Hence, A differs from A only by a coefficient. 

As a candidate for the index of biaxiality we may introduce 

«= {(Ai- A 2 )(A 2 - A 3 )(A 3 -A!)} 2 . (41) 

The index k is non-negative. It is obvious that the tensor is biaxial if and only if k =/= 0. The index k 
is called the discriminant in the context of theory of algebraic equations. We explain this point briefly. 
The eigenvalues of the matrix t = (tij) are roots of the cubic equation 

det(xl - t) = x 3 + ax 2 + bx + c = 0. (42) 

The coefficients and roots are related as 

0= -(Aj +A2 + A3) = -Tit, (43) 
b = A X A 2 + A 2 A 3 + A3A1 = ~(A X + A 2 + A3) 2 + + A 2 ) = i(Trt) 2 - ^Tr(t 2 ), (44) 

c = -AiA 2 A 3 = -dett. (45) 

k = is a necessary and sufficient condition for existence of a multiple root. It is known that the 
discriminant is expressed in terms of the coefficients as 

K = a 2 b 2 + I8abc - 46 3 - 4a 3 c - 27c 2 . (46) 

Similarly, the degree of anisotropy (|3"5|) is expressed in terms of the matrix t as 

A = (A 1 -A 2 ) 2 + (A 2 -A 3 ) 2 + (A 3 -A 1 ) 2 

= 2(A 2 + A 2 + A 2 .)-2(A 1 A 2 + A 2 A 3 +A 3 Ai) 

= 2Tr(t 2 ) - {(Trt) 2 - Tr(i 2 )} 

= 3Tr(t 2 ) - (Trt) 2 . (47) 

We may use A and n as indices of anisotropy and biaxiality. But, particularly, k is not convenient 
for calculation. What is worse, the index k is not useful for comparing the biaxiality of the microscopic 
tensor t with the biaxiality of the macroscopic tensor (t) . In the next section we will introduce a more 
convenient method to evaluate and classify anisotropies. 



4 Anisotropy diagram 

Here we will introduce a diagrammatic method to characterize anisotropy of a given tensor. Our 
diagram will be convenient for comparing anisotropies of various tensors. It will be shown that the 
microscopic tensor and the macroscopic tensor have a definite relation in our diagram. 

In the following any tensor is assumed to be traceless and symmetric. The eigenvalues of a tensor t 
are denoted as Ai, A 2 , A3. It is a usual convention to arrange the eigenvalues in the order A3 > Ai > A 2 . 
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Figure 1: In the anisotropy diagram a tensor t — diag(Ai, A2, A3) is represented by a point T = (ex, £2)- 



The lengths of segments (OHi, OH 2 , OH 3 ) are equal to \/|(Ai, A2, A3). 



Under the assumption of tracelessness the sum of the three eigenvalues is zero. Here we give definitions 
for classification of traceless symmetric tensors. The tensor is isotropic if t = 0. Otherwise, it is 
anisotropic. When two of the three eigenvalues coincide, the tensor is uniaxial. Moreover, when a 
uniaxial tensor satisfies dett > 0, namely, A3 > > Ai = A2, it is said that the tensor has positive 
uniaxiality. A positively uniaxial tensor has eigenvalues (Ai,A2,A3) — A(— 1,-1,2) with a positive 
coefficient A. On the other hand, when a uniaxial tensor satisfies dett < 0, namely, A3 = Ai > > 
A2, it is said that the tensor has negative uniaxiality. A negatively uniaxial tensor has eigenvalues 
(Ai, A2, A3) = A(l, —2, 1) with a positive coefficient A. When dett = and t ^ 0, namely, A 2 = —A3 
and Ai = 0, it is said that the tensor has maximal biaxiality. A maximally biaxial tensor has eigenvalues 
(Ai, A2, A3) = A(0, —1, 1) with a positive coefficient A. 

Any traceless symmetric tensor can be diagonalized and parameterized in the form 

(48) 

The coefficients (£1,62) are the same things as (05,04) in (f2H|) and they are related to the eigenvalues 
as 

Ai = — 9=£i + ^=e 2 = a/^ • t^(-1, V3) ( t 1 ) , (49) 




\i = -— — : 2 - A /|--(-l,-V3)^), (.-,01 

l^ = Vi (1 ' H2)- (51) 



or inversely 



ei = ^-X 3 , (52) 

£ 2 = ^(Ar - A 2 ) = -±=(2Ai + A 3 ) = ^(-2A 2 - A 3 ), (53) 
= A! + A2 + A 3 . (54) 
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Figure 2: A permutation of the eigenvalues A2 <-> A3 induces the transformation fa of the anisotropy 
diagram, which is a reflection about the line £\. Similarly, other permutations A3 <-» Ai and Ai <-> A2 
induce fa and fa, respectively. The permutations of the eigenvalues generate six equivalent points. 
There is a unique representative point in the fundamental domain F. 



The inner product (|17|) is used to define the norm of the tensor 



a = \\t\\ = v^M) = v/(Ai) 2 + (A 2 ) 2 + (A3) 2 = x/(ei) 2 + (£ 2 ) 2 - (55) 

Then a 2 is equal to the anisotropy index 3 A which was defined at (|39[) . For the tensor t we plot a 
point T whose Cartesian coordinate is (ei, £2) as shown in Fig. 1. Thus each point in the plane defines 
a corresponding traceless symmetric tensor. This plane diagram is called an anisotropy diagram. The 
value of a is equal to the distance between the point T and the origin O of the coordinate. 

We explain how to draw the anisotropy diagram in detail. For a given traceless tensor t one 
calculates the eigenvalues (Ai, A2,As). Next one calculates (ex, £2) using Eqs. ([52^1 . (fSU)) . Plot a point 
T = (81,62) in the Cartesian coordinate. This is the point representing the tensor. Draw three 
lines ^1,^2,^3 which run through the origin O = (0,0) in the direction (— 5, 3 )' (~!> — 2 )> (1>0)> 
respectively. Draw a line mi which runs through the point T and is perpendicular to the line l\. The 
intersection of £\ and mi is denoted as Hi. Similarly, draw lines m,2 and 7773 which run through the 
point T and are perpendicular to the line £2 and £3, respectively. The intersection of £2 and m-2 is 
denoted as H2. The intersection of £3 and 7713 is denoted as H3. Then the lengths of OHi, OH2, OH3 
are equal to Ai, y/| A2, y/| A3, respectively. In this way, we can determine the set of eigenvalues 
from the representing point, and vice versa. 

If it is not requested to arrange the eigenvalues in the order A3 > Ai > A2 and it is allowed to 
rearrange them, a point in the anisotropy diagram corresponding to the given tensor is not unique. The 
operation exchanging A2 «-> A3 induces a transformation of the coordinate of the anisotropy diagram 



as 



3 . /3 / 1 1 A 1 VS 

1 ^ . ^ 1 /-3 1 \ a/3 1 

71 (2Al + A2) = T2 { Te £1 + T2 £2 J = ~ £1 + 2 £2 ' 



which can be summarized as 



zi) Vi) 2 U\/3 1 y V £ 2 



(56) 
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The permutation A3 <-> Ai induces a transformation 



1 



£ i = V o A i = V o 75 £ i + -S £ 2 = ~ r £1 + £ 2, 



V3 



7i ( - 2A2 



V6 

1/3 1 A \/3 

Al) = 7lU ei + 71 £ V = ~ e 



namely, 



1/-1 V3 
2 l\/3 1 



Another permutation Ai <-» A2 induces a transformation 



2 £2 ' 



(57) 



»3 : 



1 

-1 



(58) 



A point T is moved to the point 0iT by the mapping 0i. Furthermore, it can be moved to the point 
020iT by 02, and so on. In the anisotropy diagram Fig. 2, the transformations 0i, 02, 03 are reflections 
with respect to the lines £1,^2,^3, respectively. The set of transformations {01,02,03} generates the 
third permutation group 63, which has 3! — 6 elements. Under the actions of S3 a generic point T 
in the anisotropy diagram leaves six points on its trajectory. These trajectory points Ti, • • • , Tq are 
equivalent as a representative of the tensor t. If we impose the condition A3 > Ai, Eqs. (|4"9")) and (fSTj) 
imply V3e\ > £2- Moreover, if we impose the condition Ai > A2, Eq. (f53f implies £2 > 0. Hence, if the 
eigenvalues are arranged to satisfy the conventional ordering A3 > Ai > A2, a unique representative 
point is chosen in the domain 

F = {(£ 1 ,£ 2 )|0<£ 2 < V3ei}, (59) 

which we call the fundamental domain of the anisotropy diagram. 

We can use the radius a and an angle 7 to parameterize the coordinate of the anisotropy diagram 

as 

£1= a cos 7, £2 = a sin 7. (60) 

By substituting these variables into (|4"9")) - (15ip and PT]) , we obtain an expression for the index of 
biaxiality 

I r . .,.-21 . . 

(61) 

2 



1 

K = — 

2 



£ 2 {3(£x) 2 - (£ 2 ) 2 }] 2 - \o? sin 2 7 (3 - 4 sin 2 7 ) 2 



If the angle 7 is varied, k takes the maximum value ^a 6 when sin 2 7 = j or 1. Hence the biaxiality 
becomes the maximum at 7 = f , § : if j 7f 5 ^f" , ■ On the other hand, K takes the minimum value 
when sin 2 7 = or |. Hence the biaxiality vanishes at 7 = 0, ^, ^,7r, 4f , ^p. In the anisotropy 
diagram Fig. 3, we introduce a family of half lines 



' u +x 


= {(ei 


£2) 


= ir(-l,v / 3)|r>0}, 


U- x 


= {fa 


£ 2) 


= -ir(-l,V3)|r>0}, 


U+y 




£ 2) 


= ir(-l,- v / 3)|r>0}, 


U-v 


= {{ei 




= -|r(-l,-\/3)|r>0}, 


u +z 


= {(ei 




= r(l,0)|r>0}, 




= {(ei 




= -r(l,0)|r > 0}. 



(62) 



These lines divide the anisotropy diagram into six domains. Tensors which belong to U± x , U± y , U± z 
are 

(2 \ /-l \ /-l 0\ 

(63) 
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Figure 3: The half lines {Uj} are uniaxial lines while the half lines {Bk} are maximally biaxial lines. 
The degree of the anisotropy of the tensor t is a = the length of the segment OT. The degree of the 
biaxiality is (3 = — 6. 



respectively. The half lines {Uj} are called the uniaxial li 
of half lines 



( R 

+^)— y 

B- X .+y — {{Si, Si 



k B- z ,+ x — {(ei, £2 
Then tensors which belong to B± x ^ y , B± y , TZ , B± ZtTX are 



lines. Similarly, we introduce another family 
= {(ei,e 2 ) = r(0,l)|r>0}, 
r(0, 1) \r > 0}, 
B +y ^- z = {{e u e 2 ) = |r(->/3, -1) | r > 0}, 
B- y , +z = {(ei,e 2 ) = -|r(-V5, -1) 1 r > 0}, 
r(V3,-l)|r>0}, 
■r(V3,-l)|r>0}. 



(64) 



B+x,-x = {(^ljEa) 



r Z 1 ° ° N 
± -7= -1 

^2 \0 ; 



B± z ^x 


are 


/0 





1 





Vo o 
































(65) 



The half lines {Bk} are called the maximally biaxial lines. 

Using the anisotropy diagram we can classify tensors and measure their degrees of anisotropy. An 
isotropic tensor is represented by the origin of the diagram. A uniaxial tensor is represented by three 
equivalent points on the uniaxial lines. A biaxial tensor is represented by six equivalent points and 
each point belongs to one of the six domains divided by the uniaxial lines. As shown in Fig. 3, a generic 
point T in the diagram is sandwiched between one of uniaxial half lines and one of maximally biaxial 
half lines, which are denoted as Uj and Bk- The degree of anisotropy is measured by the radius OT 
= a defined at |K5]) . The angle between Uj and Bk is (tt/6) rad. Let 9 be the magnitude of the angle 
formed by the half lines OT and Uj measured in radians. Then we can define the degree of biaxiality 
(3 by 

(66) 



7T 



Then (3 takes its value in the range < < 1. Bates and Luckhurst [23] gave another definition of 
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biaxiality index 77, which they called the relative biaxiality, 



V = — \ (67) 



A_ 

for A3 > > Ai > A 2 . It is equal to 

?? = V3^=V3tang= (68) 
ei tan(7r/6) 

The index 77 is a monotonically increasing function of [3 and it also takes its value in the range < 77 < 1. 
The intuitive meaning of 77 is also clear. As another index of biaxiality we may define 

C = t£- (69) 
a" 

The index £ also takes its value in < £ < 1 since the value of n is within < n < ^a 6 as discussed 
at (EI]). 

Here we need to mention that diagrams which are similar to our anisotropic diagram can be found 
in literatures. Kralj, Virga, and Zumer j!8j introduced a diagram which is equivalent to the anisotropic 
diagram. A new point of our study is that we use the diagram as a tool for comparing various tensorial 
quantities and for measuring the degrees of anisotropies. Another new point is that we establish a 
relation between the microscopic quantity and the macroscopic quantity in the anisotropic diagram as 
will be discussed in the next section. For a comparison with Kralj 's parameterization, Eq. (9) in their 
paper [15] . we write the eigenvalues (H9" l) - ([5Tjl in terms of the variables (I^U|) as 

Al = \Zf a C ° S ( 7 ~~ T) ' A2 = a C ° S ( 7 + T") ' A 3 = COS7. (70) 

5 Micro-macro relation 

In the previous section we introduced the anisotropy diagram, the degree of anisotropy a, and the 
degree of biaxiality f3. We have introduced also the geometric order parameter G which connects the 
microscopic tensorial quantity t with the macroscopic tensorial quantity (t) = Gt. In this section 
we will discuss how the degrees of anisotropy of the macroscopic quantity is related to the degrees 
of anisotropy of the microscopic quantity. We will establish a diagrammatic relation between them, 
which we call the micro-macro relation. 

The idea of the micro-macro relation was inspired by Ojima's idea, Micro-macro duality [27 . 
Micro-macro duality means bi-directional relations between the microscopic quantum world and the 
macroscopic classical world. Although our present consideration is restricted within classical physics, 
the relation between the molecular quantity and the macroscopic quantity can be regarded as an 
example of Micro-macro duality. 

Let us confirm notation to be used below. The orientation of a molecule in liquid crystal is described 
by a three-dimensional rotation matrix R £ SO(3). The molecules have various orientations and their 
statistical distribution is described by the probability density p(R). Assume that each molecule has a 
physical quantity which is represented by a traceless symmetric tensor t. When a molecule is turned 
by a matrix R, the tensor is transformed to (R <B> R)t = RtR T . Then the average 

(t} = j '{R® R)tp(R)dR = J {RtR T )p{R)dR = Gt (71) 

is the quantity observed at the macroscopic scale. The tensor R t R T has the same set of eigenvalues 
as t for any rotation matrix R. However, in general, the principal axes of R! t R' T do not coincide with 
those of RtR T for different matrices R and R '. In other words, the matrices {RtR T } defined with 
various R G 50(3) are not simultaneously diagonalizable. Hence, it seems nontrivial to find a general 
relation which holds between the eigenvalues of t and those of (t) . This is actually what we found and 
is called the micro-macro relation. 
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Figure 4: The points {Ti, • • • , Tg} are equivalent points representing a tensor t. The point T represents 
the average tensor (£). The point T is always in the polygon D whose vertices are {Ti, • • • , Tg}. (a) 
T is uniaxial, (b) T is generic, (c) T is maximally biaxial. 



Let {Ti, • • • , Tg} be the set of equivalent points in the anisotropy diagram which corresponds to 
the tensor t. The degrees of anisotropy and biaxiality of t are written as a(t) and f3(t), respectively. 
Let T be a point in the anisotropy diagram Fig. 4 which corresponds to (t). Let D be a convex polygon 
which has the points {Ti, • ■ ■ , Tg} as its vertices. Both the perimeter and the inner domain of D are 
included in D. Let P be an arbitrary point of D. Then the following two theorems hold. 
Theorem 1. For any probability distribution p(R), the average point T is in the polygon D. 
Theorem 2. For any point P in D, there exists a probability distribution p(R) such that the average 
point T coincides with the point P. 

Before showing proofs of these theorems, we introduce three kinds of indices which enable us to 
compare anisotropies of microscopic and macroscopic quantities. We call the ratio of the degrees of 
anisotropy 

Xi = —ttt (72) 
a(t) 

the strength of realization of anisotropy. Theorem 1 implies that the length of the segment OT cannot 
be longer than the length of OT. By definition, OT = a((t)) and OT = a(t). Hence, their ratio %i 
is in the range < Xi < 1- I n particular, for the pair of the microscopic singular tensor t a and the 
macroscopic singular tensor u a which satisfies Eq. (|3 1 [) , the strength of realization of anisotropy Xi is 
equal to the degree of order a a . 

We call the ratio of the degrees of biaxiality 

m)) (n , 
X2 = lm (73) 

the strength of realization of biaxiality. The value of xi can be larger than 1. In such a case it is 
said that biaxiality is enhanced in the macroscopic phase. It can happen that /?((£)) 7^ even when 
/3(f) = 0. In such a case we formally write xi = 00 an d say that biaxiality is generated from uniaxial 
molecules. Conversely, it also can happen that f3((t)) = although f3(t) ^ 0. In such a case we have 
Xi = and say that biaxiality is lost or uniaxiality is realized from biaxial molecules in the macroscopic 
phase. 

We define the third index 

X3 =sgn(det(t) debt), (74) 

and call it the relative signature. Here sgn (x) denotes the signature of x and it takes its value in 
{1,0,-1}. When X3 = 1, we say that the macroscopic phase is positively oriented or it has prolate 
order. When X3 = — L we say that the macroscopic phase is negatively oriented or it has oblate order. 
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We would like to show another theorem. Theorem 1 is a corollary of this theorem: 
Theorem 3. Let A max and A m i„ be the maximum and the minimum, respectively, among the eigen- 
values {Ai, A2, A3} oft. Let {/Xi, /Z2, ^3} be the eigenvalues of (t) . Then, it holds that 

Amin < Mr 5: A max , T = 1, 2, 3. (75) 

Proof of theorem 3: Let v r be a normalized eigenvector satisfying (t)v r = fi r v r . This equation is 
to be read as a multiplication of the matrix (t) on the vector v r as 

3 

2J (tij)Mj = Vr{v r )i- (76) 

We write an inner product of vectors u and w as (v.w) — ^2i=i w i ^i- Then we have 

= (v r ,(t)v r ) 

'v r ,J (Rt R T ) p(R) dRv r ^ 

(R T v r ,tR T v r )p{R) dR 

= J Trjt (i? T v r i? T v r ) | p(i?) dR 

= Trjf y (i? T v r ® i? T v r ) p(i?) di?} 

= Tr(tp r ). (77) 
In the last line we introduced the three-dimensional matrix p r which is defined by 

p r = J {R T v r ® R T v r )p{R) dR = J R T (v r ® v r )Rp(R) dR. (78) 

The matrix p r is symmetric, non- negative and satisfies Tt p r = 1. On the other hand, let 

t = Aiir + A 2 n 2 + A 3 n 3 (79) 

be the spectral decomposition oft. The three-dimensional matrices {n 9 } satisfy (n g ) T = II g , n^IIs = 
S qs U s , Y/ q =i n 9 = I, TrlL, = 1. Substituting this into dT7J> we obtain 

3 3 

p r = Tr(p r t) = 53 Tr(p,n 9 ) A 9 = ^ P r , A g (80) 

g=l g=l 

where P rg = Tr(p r IIg) are non-negative real numbers and satisfy frg = 1< Hence, 

3 3 

Mt" ^ ^ ^7~g A^ ^ ^ ^7~g A max A max , ($1) 
q=l q=l 
3 3 

/l r ^ ^ ^rg A^ ^ ^ ^ -frg A m j n A m i n . ($^) 
g=l g=l 

This ends the proof of theorem 3. It is also interesting to note that J2r=i Pr = I an d therefore 
E'=i^g = l- 

Proof of theorem 1: From the construction of the anisotropy diagram it is obvious that the set of 
points satisfying the inequality ((75]) is the polygon D. See the Fig. 5. This fact can be also verified 
via explicit calculations. The coordinate (£1,62) of the point T is defined by 

/3 1 
2 M3 ' ^2 = ^(^1-^2)- (83) 
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Figure 5: The intersection of the shaded bands, Zx, Z 2 , Z3 defines the polygon D. The average point 
is restricted in D. 



They satisfy the equation similar to (|4"9"l) - (fB"Tj) . 



1 . 1 _ 12 1 , . ~, ie x 



(-1,n/3) [V- , (84) 



^#-v^ (1 < 0) (S)- (86) 

Hence the sets of points restricted by the inequality (f?5|) . 

Z\ = {(ex,e 2 ) I A min < --^= ex + -j= e 2 < A max }, (87) 

Z 2 = {(El, £2) I A m in < — — ^= £i — — /= £2 < A max }, (88) 

Z 3 = {(ei, £2) I A min < ^1 ei < A max }, (89) 

are drawn as three shaded bands in Fig. 5. Their intersection Zx H Z 2 n Z3 is nothing but the polygon 
D. This observation proves theorem 1. 

Proof of theorem 2: First, let us note the following simple fact. Suppose that two traceless symmet- 
ric tensors to an d ti are simultaneously diagonalizable. In other words, they have common principal 
axes. Let To and Ti be their representing points in the anisotropy diagram. Then the weighted sum 

t w = (1 - w)t + wtx (90) 

with a real number w (0 < w < 1) is also a traceless symmetric tensor and diagonalizable simultaneously 
with to and t\. It is easily verified that the point representing t w in the anisotropy diagram divides 
the segment T Ti in the ratio w : (1 — w). 

Second, let us note that any point P of a convex polygon can be expressed as a weighted sum of 
the vertices of the polygon. We can chose real numbers Wx, • ■ ■ ,wq such that 

6 6 

t P = s jTw l t l , J2 w i = 1 > 0< Wl <l (91) 
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Figure 6: The point Q divides the edge T3T4 into the ratio 1 : 2. The point P divides the segment 
TiQ into the ratio 3 : 1. Thus, Q = §T 3 + ±T 4 , P = ±Ti + §Q = ±Ti + ±T 3 + ±T 4 . Every point in 
the polygon D can be expressed as a weight sum of the vertices. 



For example, the point Q in Fig. 6 is given by 

*Q = ~*3 + (92) 

and the point P is given by 

13 111 
fp = -h + ~t Q = ~t x + -t 3 + -U. (93) 

For a given point P the set of weights (wi, ■ ■ ■ ,We) is not unique but uniqueness is not necessary. 

Third, remember that the vertices of the polygon D in the anisotropy diagram are related to each 
other by reflections fa, fa, fa. Note that these reflections in the diagram can be generated by rotations 
in the real space. If we define 

/l \ / l\ A) -1 0\ 

#1=00-1, # 2 = 10, #3=1 , (94) 
\0 1 / \-l 0/ \0 1/ 

then the reflection mapping fa introduced in Sec. 2] is equivalent to the rotation fat = Ki tKj , Hence, 
for the vertices t\(= t), ti, • ■ ■ , *6 of the polygon there exists a set of rotation matrices Ri, R2, ■ ■ ■ , Re 
such that ti = Rit Rj . 

Finally, combining the above arguments we obtain 

6 6 

t v = ^w i t i = Y J w i {RitRj), (95) 

i=l i=l 

which should be compared with Eq. (ffTj) . This means that the set of weights [w\,--- , We) is a 
probability distribution which yields the point P as the average. This proves theorem 2. 
The collection of theorem 1,2, and 3 is called the micro-macro relation. 



6 Examples 

In this section we will demonstrate calculations of the macroscopic tensors by assuming simple prob- 
ability distributions. In the first example we will show that uniaxial molecules can generate a biaxial 
order at the macroscopic scale. In the second example we will show that uniaxial molecules can ex- 
hibit a negatively oriented uniaxial phase. The third example has a continuous probability distribution 
and will exhibit the same result as the second one. In the fourth example we will show that biaxial 
molecules can exhibit a uniaxial order. 

In the first example we assume that the molecule has a physical quantity 

(96) 
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Figure 7: Examples of calculation of the average tensor T from the microscopic tensor T. (a) The 
uniaxial tensor T can generate the biaxial TV It can also generate the uniaxial T2 with the inverted 
signature, (b) The biaxial tensor U can generate the uniaxial U. 




which has positive uniaxiality on the z-axis. The anisotropy diagram for this is shown in Fig. 7 (a). 
The coordinate of the representing point T is (£1,22) = (1,0). Assume that 2/3 of molecules are 
aligned in the ^-direction and 1/3 of molecules turn into the n-direction. Then the average of the 
tensorial quantity is 

2 1 2 1 (- 1 °\ 1 1 ( 2 

(97) 

The coordinate of the point Ti corresponding to (t)i is (ei, £2) = |(1, and it lies on the half line 
B-y,+z- The maximal biaxiality is realized in this case. The strength of realization of anisotropy is 

Xi = ^ = /=- (98) 
Xl a(t) V6 V3' ( j 

The degrees of biaxiality are f}(t) = and /3((t)i) = 1. Hence the strength of realization of biaxiality 
is X2 = 00. 

Next, assume that 1/2 of molecules turn in the ^-direction and 1/2 of molecules turn into the 
y-direction. Then the average is 

I I 1 1 I 2 °\ I I f- 1 

(t) 2 =-(K 2 tKZ) + -(K 1 tK?) = - ■— -1 +-._ 2 



, n n -i / 

2 76 




-1, 



(99) 



The coordinate of the point T2 in Fig. 7 (a) corresponding to (2)2 is (£1,62) = (— 1,0). The strength 
of realization of anisotropy is \i — \- The degrees of biaxiality remains f3((t)2) — 0. The relative 
signature is X3 = — 1 m this case. Thus the macroscopic phase has oblate order. 

In the above two examples, the probability distribution p{R) that we assumed had pointwise sup- 
port, namely, the integral in Eq. (|71[) was replaced by summation. Here we show an example which 
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has a continuous probability distribution. We define a rotation matrix 




K :l [o) ■■ ( sin.. <•«,*.> j (100) 

which is parameterized by an angle 0. Assume that the molecules are turned as 1 1— > Ks(c/))K2 1 K^K^tf)) 7 
with a probability distribution which is uniform with respect to the variable <p. Then the average be- 
comes 

r 27T 



(K3(4>)K 2 t K2 K 3 ((j)) ) d<j) 



'0 




2 cos 2 4> — sin 2 </> 3 cos 4> sin (/) 
3 cos <p sin <^> 2 sin 2 — cos 2 < 





(101) 



This result is the same as (|99|) . 

In the fourth example we assume that the molecule has a physical quantity 




u= —j= -1 . (102) 



The representing point U is shown in Fig. 7 (b). Its coordinate is (£1,^2) = ^(\/3, 1)- It lies on the 
half line B- Vy+Z and has the maximum biaxiality (3{u) = 1. Assume that 1/2 of molecules are aligned 
in the same orientation and 1/2 of molecules are turned about the z-axis by the right angle. Then the 
average becomes 




(103) 

which has positive uniaxiality on the z-axis. In this case biaxiality is lost in the macroscopic phase. The 
coordinate of the representing point U is (e"i, £2) = "2 (1 ; 0). The strength of realization of anisotropy 

is Xi = # 



7 Order parameters for Z^ft-symmetric molecules 

Here we apply the method of the geometric order parameter to platelet molecules. By platelet molecules 
we mean molecules which possess the point group D211 as its symmetry. The shape of a platelet molecule 
is invariant under reflections on the yz-plane, xz-plane and xy-plane. The T>2h symmetry is a group 
generated by 

-1 0\ /l 0\ /l \ 

10 , 0-10, 01 . (104) 

1/ \0 1/ \0 -1/ 

A physical quantity t of the platelet molecule must be invariant under the actions of D2h- Namely, 
it is required that RtR T = t for any R £ Dih- This implies that the off-diagonal elements satisfy 
— tij = tij for i ^ j, hence Uj = 0. Thus, only the diagonal elements ta can be nonzero. Moreover, 
we assume that the associated macroscopic quantity (t a b) = Si j Gabijtij is also invariant under the 
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actions of D<zh- Namely, it is assumed that R(t)R T = (t) for any R £ P>ik- Then only the diagonal 
elements (t aa ) can be nonzero. The invariance under the D^h transformations requires that 



3 3 
^ Rac Rbd Gcdij — G a bij , ^ G a bij RikRji = G a bki for (Rij) S D2h- (105) 

c,d—l iji=l 

These requirements for the geometric order parameter G are equivalent to 

- Gabij = Gabi-j for a^b, ~G ablj = G abio for i ^ j (106) 

Therefore, 

Gabi] =0 for a ^ b, G abi] = for i ^ j. (107) 

Hence, elements which can be nonzero are G aa n with a,i = 1,2,3. In the following we abbreviate it 
as G aa ii = G a> i- The nine components {G a -i} are imposed the traceless condition 

3 3 

^G a;J =0, ^G a;i = 0. (108) 

a— 1 i— 1 

Hence, only four components among {G a -i} are independent. 

Dummur and Toriyama [HI [H] dehned four parameters S,D,P,C as 

(109) 

^(G 3;1 - G 3;2 ), (110) 
= ^(G 1;3 -G 2;3 ), (111) 



5 = 


c-ZZ 
°zz — 


2 G3;3 


L> = 


qZZ _ 


qZZ _ 
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P = 


qXX 

°zz 


qYY 

°zz ' 


G = 


( oXX 


qYY 

^xx 



yy yy 



;2 + G 2;2 ). (112) 



These are almost equal to the parameters S, T, U, V which Straley [5] introduced. Only the differences 
between (S, D,P,C) and (S, T, U, V) are multiplicative factors as explained in the reference [24]. Here 
Sff are elements of the de Gennes ordering matrix The indices a,b = X,Y,Z specify axes of 
the laboratory frame while i,j = x,y,z specify axes of the molecular frame. Physical meanings of 
these parameters are explained as follows. The parameter S is an index to measure how strongly the 
uniaxiality of the molecule manifests itself as the uniaxiality of the macroscopic phase. D measures 
how strongly the biaxiality of the molecule manifests itself as the uniaxiality of the macroscopic phase. 
P is an index of the strength with which the molecular uniaxiality generates the macroscopic biaxiality. 
C represents how strongly the molecular biaxiality generates the macroscopic biaxiality. 

The parameters (5, D, P, C) can be written in various forms. In the present case the microscopic 
quantity t and the macroscopic quantity (t) can be diagonalized as t = diag(Ai, A 2 , A 3 ) and (t) = 
diag(/ii, /i 2 , /i 3 ). Then the relation (t) = Gt is written as 

1:3\ /Ai\ 

(113) 

The set of equations (|108p . p09p - (|112|) is solved for the geometric order parameters as 

'Gi ; i Gi ;2 Gi ;3 \ 1 /S-D-P + C S + D-P-C -2S ■ 
Gaji G 2 ; 2 G 2 ; 3 =- [S-D + P-C S + D + P + C -2S-2P]. (114) 



Gl;l 


Gl ;2 


Gi ;3 \ 


fAi 


G 2 ;l 


G 2;2 


G 2 ;3 j 


A 2 


G3;l 


G3 ;2 


G 3 ; 3 / 


^A 3 




v G 3; i G 3;2 G 3;3 / V -2S + 2D -2S-2D 
Eqs. (jTT5f and (TTT4"|) yield 

M3 = \{~S+ D)M + \(S- D)X 2 + ~S\ 3 = 1d(Ai - A 2 ) + SA 3 , (115) 

Mi - M2 = \(-P + C)Ai + \(-P- G)A 2 + ^PA 3 = i(7(A x - A 2 ) + PA 3 , (116) 
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where we used Ai + A2 + A3 = 0. They can be put in the form 



This equation is consistent with the interpretation of the parameters {S, D, P. C) explained above. It 
can be expressed in terms of the anisotropy coordinates ([52l . (|53|) . (|83|) as 

(118) 



£2) UP \C J \e 2 ) {d i5 6 J \e 2 



On the other hand, the elements of the reduced ordering matrix G are calculated from the definitions 
(HU), (HU) with the help of (fTUS)) as 

1 3 

G55 = —{Gl-l + Gi : 2 — 2G*i : 3 + G^-A + G 2 ;2 _ 2G2 ; 3 — 2G3 ; i — 2G3 ; 2 + 4G3 ; 3) = — Gs-3 = S, 

b 2 

1 1 

G54 = — /={" Gi-i + Gi ; 2 — G^-i + G2 ; 2 + 2G3 ; i — 2G3 ; 2) = —^-{Gs-i — G?,^) = 

1 i 

G45 = — f= i~ Gl;l — Gl-2 + 2Gl;3 + G2;l + G2;2 — 2G2;3) = -^-(Gl;3 — G2;3) = ^J^P' 

G44 = -(Gi ; i — Gi ; 2 — G2 ; i + G2 ; 2) = —C 
This result is consistent with (|1 18|) . 



8 Landau-de Gennes free energy 

The Landau-de Gennes free energy is a standard tool for analysis of phase structures of liquid crystals. 
The Landau-de Gennes free energy is a polynomial function J^(a) of a collection of macroscopic 
quantities, which is denoted as a. The quantities a play the role of order parameters, too. The free 
energy should be invariant under spatial rotations of the variables. This requirement is symbolically 
written as J^(Ra) — ,^{a). The coefficients ci, C2, • • • in the polynomial J^"(a) = cia + c-2a 2 + ■ ■ ■ may 
depend on various external physical parameters like temperature or density. It is required that in an 
equilibrium state the free energy takes its minimum value. Thus the values of the order parameters a 
are determined as the minimizer of the free energy. Then the symmetry of the equilibrium phase is 
determined by the values of the order parameters, which are also functions of the external parameters. 
This is a usual routine to analyze the phase structure using the Landau-de Gennes free energy. 

In this section we will explain a general prescription to formulate the Landau-de Gennes free energy 
of arbitrary shape molecules. Here we mainly consider nematic phases, which are translationally 
invariant. A possible generalization for smectic phases will be discussed briefly. Later we will apply 
our prescription to the Z^-symmetric molecules. 

Assume that each molecule has microscopic quantities ao, bo, • • • , which are symmetric tensors. 
They can be the dielectric susceptibility tensor or the magnetic susceptibility tensor of a molecule. 
We can assume that they are traceless. If cti is not traceless, we can take the traceless component 
by subtracting its trace to define — a\ ~ iTr(ai). At a macroscopic scale we measure physical 
quantities, a, b, ■ ■ ■ , which are ensemble averages of microscopic quantities. The geometric order 
parameter G relates them as Gao = a. The tensor a is transformed under a rotation R 6 5*0(3) as 
a 1— > RaR T . A polynomial function I(a, b, ■ ■ ■ ) which satisfies 

I(RaR T : RbR T , ■■■) = I(a, b, ■ ■ ■ ) (119) 

for arbitrary R £ 5*0(3) is called an invariant polynomial. Then the Landau-de Gennes free energy is 
defined as a function of G, 

&{G) = I(Ga , Gb Q , ■ ■ ■). (120) 
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The invariant polynomials up to the fourth order are listed as 



h 


= Tr(a&) = Tr(6a), 


(121) 


h 


= Tr(abc) = Tr(acb), 


(122) 


h 


= Tr (abed) = Tr (deba), 


(123) 


h 


= Tr (abdc) = Tr (cdba), 


(124) 


h 


= Tr (acbd) = Tr (dbca), 


(125) 


h 


= Tr(ab)Tr(cd), 


(126) 


h 


= Tr (ac) Tr(6d), 


(127) 


h 


= Tr(ad)Tr(bc). 


(128) 



The equal signs in the above equations hold since a T = a and 

Tr(abc) = Tr(abc) T = Tr(c T b T a T ) = Tr{cba) = Tr(ac6). (129) 

Other polynomials like Tr(a)Tr(6), Tr(a)Tr(foc), Tr(a)Tr(6) Tr(cd) become zero because a, b, • • • are 
traceless. Using the representation theory of the rotation group, we can prove the following theorem: 
Theorem 4. Among the polynomials of traceless symmetric tensors, there is only one linearly inde- 
pendent invariant of the second order, that is I\ . There is one linearly independent invariant of the 
third order, that is 1 2 . There are five linearly independent invariants of the fourth order. The six 
invariants {I3, ■ ■ ■ ,1%} always satisfy 

h+h + h-\(h+I 7 + h) = 0- (130) 

Hence, only five among {I^,--- ,1%} are linearly independent. 

A proof of this theorem is given in the appendix. If there are four independent microscopic quantities 
Gto, bo j c 0: do, the Landau-de Gennes free energy up to the fourth order is constructed by substituting 
a = Gao, b = Gbo, c = Gcq, d = Gdo with possible repetitions into (I12ip - (|127j) and by making their 
linear combinations as 

&(G) = c 1 I 1 (Gao,Gao)+C2li(Gao,Gbo)+c 3 Ii{Gao,Gco) + c i Ii(Gao,Gdo) 
+c 5 h(Gb , Gb ) + c 6 h(Gb , Gc ) + c 7 h{Gb Q , Gd Q ) 
+c B I 1 (Gc , Gc ) + cqIx(Gco, Gd ) + ci ii(Gd , Gd ) 
+c 11 I 2 {Ga () , Ga , Ga ) + ci2/ 2 (Ga , Ga , G6 ) + ci 3 / 2 (Gao, Ga , Gc ) 
+c 1A I 2 (Ga Q ,Ga Q , Gd Q ) + c 15 I 2 (Ga , Gb Q , Gb ) + ■■■ . (131) 

The coefficients ci, c 2 , ■ ■ ■ may depend on temperature or density of the liquid crystal. The values of 
the order parameters, G a Mj or G M1 ,, are determined as the solution of the minimization problem of the 
free energy J^"(G). However, it can happen that the solution G a bij take physically unrealizable values. 
The degrees of order a a , which were define at (|3ip . must take their values in the range < a a < 1 
to be physically realizable. If o~ a is larger than unity, the values of the geometric order parameters 
determined by the Landau-de Gennes free energy model should be regarded as an unphysical wrong 
solution. 

When an external electric field or magnetic field is applied, the rotational invariance is broken and 
hence the free energy can have extra terms. If the molecule has an dielectric susceptibility tensor 
s = (sij) and if an electric field E is applied, the free energy has an additional term 

3 

J? (G) = E ■ Gs ■ E = E aE b Gabij . (132) 

a,b,i,j=X 
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On the other hand, if the molecule has an electric quadrupole moment q = (qij) and if an inhomoge- 
neous electric field E(x) is applied, the free energy gets an additional term 

BP 

& q (G)=Tr(\7E-Gq) = £ G a bij qij- (133) 

a,b,i,j — l 

In most of our discussions we are treating only nematic phases. Here we briefly discuss other phases 
which are not translationally invariant. In smectic or cholesteric phases, the tensorial quantity t(x) = 
G(x)to can depend on the space coordinate x = (x\, x 2 , £3) and the free energy of a continuum model 
has an extra term which is expressed as a spacial integral 

3 



j3 x 



& \G\ = [ ( k + k d^atA d 3 

k J , ,\ 1 dx a dx a 2 dx a dx b J 

a,o,c— 1 x ' 

= / E ( fci ^f ^Sf +fc2 ^Sf ^)*°«*o#**- ( 134 ) 

Then the free energy becomes a functional of G a bij(x). 

Let us apply our general scheme to a system which consists of Z^/i-symmetric molecules. The 
point group D 2 h is generated by the set of transformations (|104|) . Any microscopic quantity a Q of a 
^2/i-symmetric molecule should satisfy Ra$R T — ao for R € -D2/i- There are only two independent 
quantities satisfying this condition, 

1 /- 1 °\ /3 A 0\ 

a = _ -10 , b = J- -1 . (135) 
^ 6 \ 2/ V 2 \0 0/ 

The quantity ao is assigned a coordinate (ei,£2) = (1,0) in the anisotropy diagram. According to 
Eq. (|118|) . the coordinate of the corresponding macroscopic quantity a — Ga is (£1,62) = (S, ^j-P)- 
Similarly, the microscopic quantity 60 is (ei,£a) — (0,v3) and the macroscopic quantity b = Gbo is 
[Si, £2) = (D,^C). They can also be expressed as 

, f-S + P 

a = —\ -S-P 01, b=-pl -D-C I. (136) 



V6 








Invariant polynomials formed with these tensors are 



Tr (ad) = S 2 + \p 2 , (137) 

Tt{ab) = SD+-PC, (138) 

Tr(bb) = D 2 + ^C 2 , (139) 

Tr(aaa) = ^(S* 3 - SP 2 ), (140) 
V6 

Tr(aab) = -^-={ZS 2 D - P 2 D - 2SPC), (141) 
3v6 

Tr(abb) = — L(3SX> 2 - S*C 2 - 2PDC), (142) 
3V6 

Tr(bbb) = -^(L> 3 - DC 2 ). (143) 
V6 



The relation (|130| implies 



Tr(aaaa) - -Tr(aa) Tr(aa) = 0, (144) 
Tr(aaab) - ^Tr(aa) Tr(ab) = 0, (145) 
2Tr(aabb) + Tr(abab) - Tr(ab) Tr(ab) - ^Tr(aa) Tr(b6) = 0. (146) 
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Furthermore, since product of diagonal tensors a, b is commutative as ba = ab, it holds that 



Ti(aabb) = Tr(abab). (147) 
Hence there are only six independent invariants of the fourth order 

Tr(aa)Tr(aa), Tr(aa)Tr(ab), Tr(aa)Tr(b6), 

Tr(ab)Tr(ab), Tr(a6)Tr(66), Tr(bb)Tr(66). (148) 

The Landau-de Gennes free energy of Z^h-symmetric molecules is constructed as a linear combination 
of these invariants, 

&{S, P, D, C) = ciTr(aa) + c 2 Tr(a6) + c 3 Tr(&6) 

+c 4 Tr(aaa) + c 5 Tr(aab) + c 6 Tr{abb) + c 7 Tr(bbb) 
+c 8 Tr(aa)Tr(aa) + c 9 Tr(aa)Tr(a6) + c 10 Tr(aa)Tr(bb) 
+c u Tr(a6)Tr(a6) + c 12 Tr(a6)Tr(66) + ci 3 Tr(66)Tr(bb) 

= Cl (S 2 + ip 2 ) + c 2 (SD + ipC) + c 3 (£ 2 + ^C 2 ) 
+-^(5 3 - SP 2 ) + -^=(3S' 2 Z3 - P 2 D - 2SPC) 

2 or>r\r^\ i C7 I n3 n/^2\ 



- SC Z - 2PDC) + -^(D d - DC 1 ) 



3V6 " V6 

+c 8 (S 2 + ip 2 ) 2 + C9 (S 2 + ^P 2 )(SD + ipC) + c 10 (S 2 + \p 2 )(D 2 + ^C 2 ) 

+c n (SD + ipC) 2 + c 12 (SD + ^PC)(D 2 + ^C 2 ) + c 13 (D 2 + ^C 2 ) 2 . (149) 

up to the fourth order. This contains more terms than the free energy formulated by Allender et 
al. [TU1 fTTj even if only terms lower than fifth order are compared. Calculation of higher order terms is 
cumbersome but feasible. The values of the order parameters S, P, D, C are determined as a solution 
of the minimization problem of the free energy ^(S, P, D, C). It should be checked whether these 
values are physically realizable or not. We can calculate singular values a a (a = 1, 2) of the reduced 
ordering matrix G, which was defined at (|118p . 



fG 55 


G5i\ _ 






\Ga5 


G44J 




\G ) 



The singular values a a should be in the range |<r a | < 1. If they are not in this range, the values of 
S, P, D, C are physically unrealizable. In such a case, the coefficients ci, C2, • ■ • in the polynomial & 
should be re-adjusted. 

It should be noted that the solution of the minimization problem of &{S, P, D, C) is not unique. As 
a remnant of the rotational symmetry of the free energy, the solutions have the permutation symmetry 
S3, which is generated by the rotation matrices {K±, K2, K3} given at (|9~4"|) . or by the transformations of 
the anisotropy diagram 4>2, 4>3} given at (|56[) - (|55|) . If (S, P, D, C) is a solution of the minimization 
problem, (S", P', D', C') which is defined by 

= 0iG = 





-S-P — ^D-^C\ 

(151) 

V3 A ^^3 V 



_ 71 7S 

-v^S + ^p -D + \C 



is also a solution. In this way we obtain a complete set of equivalent solutions, {G, <j>iG, foG, </> 3 G, 
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020iG, (frifoG}, which are written as 



(S',P',D',C') = (S,P,D,C), 

-(-S - P, -35 + P,-D- C, -3D + C), 

-(-S + P,3S + P, -D + C,3D + C), 
(S,-P,D,-C), 

-(-S + P, -35 - P, —D + C, -3D - C), 
-{-S-P,3S-P,-D-C,3D-C). (152) 

9 Conclusion 

In the introduction of this paper we pointed out that the conventional method using tensorial order 
parameters ([1]), ([3]), ([5]) for characterizing biaxial nematic phases becomes ambiguous when it is applied 
to a system of asymmetric molecules. Since the conventional tensorial order parameters depend on the 
choice of a reference frame fixed on the molecule and an asymmetric molecule does not have preferable 
axes, the order parameters are not defined uniquely. What is worse, an asymmetric molecule may 
possess various tensorial physical quantities which do not have common principal axes. Although the 
ordering matrix which was originally introduced by de Gennes, is applicable to a molecule which 
has an arbitrary shape, the interpretation of the ordering matrix is difficult. Thus, we aimed to invent 
useful tools for describing and for analyzing geometric structures of biaxial nematics. 

Here we summarize the main results of this paper. Around Eq. (|15|) we argued that the ordering 
matrix is to be understood as the geometric order parameter G a uj = {QaUj) which relates the mi- 
croscopic quantity t intrinsic in a molecule to the macroscopic quantity (t) = Gt observed in a bulk 
system. The geometric order parameter was analyzed by the singular value decomposition. At Eq. pip 
it was shown that the microscopic singular tensor t a manifests itself as the macroscopic singular tensor 
u a in the nematic phase with the strength of realization a a . 

Any tensorial quantity is mapped in the anisotropy diagram. It should be noted that six or 
three equivalent points in the anisotropy diagram correspond to one tensor. As indices for evaluating 
anisotropics of tensorial quantities, we introduced the degree of anisotropy a at Eq. (|55p and the degree 
of biaxiality (3 at Eq. (|66[) . The index a is the radius in the anisotropy diagram and (3 is the angle 
measured from the uniaxial line in the anisotropy diagram. By proving theorems 1, 2 and 3 we showed 
the micro-macro relation, which tells that the point representing the macroscopic tensor always locates 
in the polygon in the diagram whose vertices are points representing the microscopic tensor. 

In Sect. [7] we applied our method to a system which consists of Z^/i-symmetric molecules. All 
the tensorial quantities of a Z?2/i-symmetric molecule have common principal axes and hence they are 
simultaneously diagonalizable. The geometric order parameter also becomes diagonal, namely, only the 
components G aa n with a, % = 1,2,3 can be nonzero. Hence it has only four independent components 
as shown in (|114p . In this case the micro-macro quantities are related as ()118p . 

In Sect. [S] we explained the general prescription to formulate the Landau-de Gennes free energy. By 
this prescription we can construct the free energy which contains all the symmetry-admissible terms 
but contains no redundant terms. We wrote down the concrete Landau-de Gennes free energy (|149p 
for the Z?2/t-symmetric molecule system. 

We would like to emphasize that we made the implication of de Gennes's ordering matrix clear 
by interpreting it as the geometric order parameter which transforms a microscopic quantity to a 
macroscopic quantity. The anisotropy diagram and the anisotropy indices which we introduced are 
systematic tool and help us understand the properties of biaxial nematics. It also should be emphasized 
that our method has no ambiguity. As an example to show the uniqueness of our procedure, we 
formulated the Landau-de Gennes free energy in the most general form without redundancy. 
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Here we mention some remaining problems. We should analyze the phase structure using the 
Landau-de Gennes free energy. This problem will be discussed in the next work. Our method can 
be applied also for analysis of molecular dynamics simulation of nematics. It can be generalized for 
treating a system which is a mixture of rod-shaped molecules and disk-shaped molecules. It can be 
generalized for treating flexible molecules although in our discussion molecules were assumed to be 
rigid. It is also interesting to study smectic phases using the continuum model which was briefly 
discussed at (|134l) . Kimura [5] studied a system of rod-shaped molecules which interact each other via 
both the short-range exclusion force and the long-range dispersion force. For an asymmetric molecule 
the principal axes of its electric quadrupole tensor may not coincide with the geometric axes which 
characterize the rigid-body repulsive force. It seems interesting to study a system which consists of 
such asymmetric molecules. In this paper we analyzed only the second-rank tensors. It is possible 
to extend our argument to include higher-rank tensorial quantities although necessary for such an 
extension is not obvious in the context of liquid crystal physics. 
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A Proof of theorem 4 

Here we prove that Eqs. (|121J) - (|128|) are a complete list of invariant polynomials up to the fourth order. 
We also prove Eq. ![T5UJ) . 



First, we count linearly independent invariant polynomials formed by products of traceless symmetric 
tensors. The set of the whole traceless symmetric tensors becomes a five-dimensional irreducible 
representation space of the rotation group. The symbol 5 or 3 denotes a five-dimensional or three- 
dimensional irreducible representation space, respectively. The one-dimensional representation space 1 
is a set of quantities which are invariant under the action of the rotation group. In other words, 1 is a 
set of scalars. According to the Clebsch-Gordan law [E], the tensor prodct space 5 (g) 5 is decomposed 
as 



In the decomposition the one-dimensional representation 1 appears once, which corresponds to I\ = 
Tr(afo) of Eq. (|121|) . Similarly, the three-fold tensor product 5 ® 5 <g> 5 is decomposed as 



h+h + h--{h + h + h) = Q. 



(153) 



5 05= 103050709. 



(154) 



5 <g> 5 <g> 5 = 



5© 



30507© 



I©3ffi5ffi7ffi9ffi 



3©5©7©9©11© 



5 709011013. 



(155) 
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In this decomposition 1 appears only once, which corresponds to I2 = Tr(abc) of Eq. (|122[) . The 
calculation of the four-fold tensor product yields 



5050505 = 1©3©5©7©9© 
3©5©7ffi 
I©3ffi5ffi7ffi9ffi 
30507090110 
50 

30507© 

1030507090 

3050709011© 

5©7©9ffillffil3ffi 

30507© 

103050709© 

3©5©7ffi9ffillffi 

507090 110 13© 

7090110130150 

1030507090 

30507090110 

507090 110 13© 

7© 9 ©11 ©13 ©15© 

9 ©11 ©13 ©15 ©17. (156) 

In this decomposition 1 appears five times. Thus there must be five linearly independent invariants 
of the fourth order and there are no more than five. The six quantities listed in Eqs. (|123p - (|128[) 
are invariants of the fourth order. By construction it is obvious that there arc no more independent 
polynomials of the fourth order. Hence the six quantities must have one nontrivial relation. From 
their symmetry, we can guess a relation of the form 

I 3 + I 4 + h + c(h 

If we substitute 



we get 1^ = 14 = I5 = 18 and = I7 = Ig = 36. Hence the coefficient c in l|157|) must be c = — |. 
This proves (|130[) . If another evidence is requested, we may substitute 



+ h + h) = 0. 


(157) 


-1 0\ 




-10, 


(158) 


2/ 






(159) 



Then we get I3 = I4 = I5 = 2 and I§ = I7 = I& = 4. This confirms that c = — h. The reader may 
calculate other cases to confirm (|153|) . 
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